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The author argues that the calculated masses of heavy tetraquarks obtained by solution of the 
spin-independent homogeneous Lippmann-Schwinger integral equation in a diquark-antidiquark pic¬ 
ture reported by M. Monemzadeh et al., Phys. Lett. B 741, 124 (2015), are incorrect. We have 
reexamined all of the published results and we believe that not only the reported tetraquark masses 
for states with zero total angular momentum are incorrect, the reported masses for states with 
non-zero total angular momentum are quite misleading, because these states cannot be predicted 
by a spin-independent formalism. 


In a recent letter by Monemzadeh et al. [H, the 
tetraquark bound state is studied as a two-body prob¬ 
lem in a diquark-antidiquark picture. The nonrelativis- 
tic s—wave bound state of a diquark-antidiquark system 
with reduced mass Pup, and the binding energy Ep, in¬ 
teracting with pair force V(r) can be described in con¬ 
figuration space by homogeneous Lippmann-Schwinger 
(LS) integral equation 
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a reduction in the calculated masses Clearly the re¬ 
ported results in Ref. [l| can not be trusted, because 
as it is shown in Table HI not only J = 0 states are not 
calculated correctly, they have reported the masses for 
non zero J states, like and AA, which cannot 

be obtained in this formalism. It is a serious challenge 
and the authors should clarify how this spin-independent 
formalism can distinguish different spin states and con¬ 
sequently the solution of spin-independent LS integral 
equation can predict the masses of tetraquarks with non 
zero total angular momentum. 


where k = \/2 Pnh I \■ The corresponding equation in 
published letter [l| (i.e. Eq. 2 and consequently Eqs. 4, 
9 and 11) is missing a factor of which should 

be a typo, because our numerical analysis demonstrates 
that missing this factor leads to completely unreasonable 
tetraquark masse^. 

In order to study tetraquark bound states, the spin- 
independent part of the diquark-antidiquark potential 
given in Ref. Q is used, with the strong coupling con¬ 
stant as of Ref. 0 
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with /3o = 11 — In/ {rif is the number of flavor quarks), 
Mo = 2.24\/A = 0.95 GeV and A = 0.413 GeV. 

In Table U we have listed our numerical results for Is 
heavy charm tetraquarks for states with total angular 
momentum equal to zero. As it is shown our calculated 
masses are all larger than the predicted results of Ref. @ 
which are obtained from solution of relativistic and spin- 
dependent LS integral equation in momentum space. Our 
results seem to be quite reasonable, because our numer¬ 
ical analysis indicates that the relativistic effects lead to 


^ To verify derivation of LS integral equation o and test our 
numerical algorithm and code, we have solved it for spin- 
independent Malfliet-Tjon Yukawa potential Q] and we have re¬ 
produced the same deuteron binding energy obtained from solu¬ 
tion of LS equation in momentum space. 


TABLE 1. Masses of charm diquark-antidiquark for Is state 
with total angular momentum J = 0 calculated from non- 
relativistic Lippmann-Schwinger equation. S and A denote 
the scalar and axial vector diquarks. 


Diquark content 


Mass [GeV] 




cqcq 



present 

Ref. ill 

Ref. 

SS 

3.885 

3.70314 

3.812 

AA 

4.013 

3.83908 

3.852 



cscs 



present 

Ref. IJJ 

Ref. 

SS 

4.117 

4.05390 

4.051 

AA 

4.250 

4.09962 

4.110 


Moreover, from numerical point of view the following 
issue should be considered. Since a regularized form of 
the diquark-antidiquark potential is used to overcome the 
singularity of the confining potential at large distances, 
the calculated masses should be independent of the reg¬ 
ularization cutoff Tc- In Ref. [H the authors have chosen 
the non-zero root of the potential as regularization cut¬ 
off and consequently the potential is fixed equal to zero 
for r > Tc- As we have shown in Table El clearly this 
cutoff is not high enough and in order to achieve the 
cutoff-independent results, converged with four signifi¬ 
cant digits, one needs to choose a cutoff at least equal 
to 5 GeV“^. Choosing a regularization cutoff equal to 
non zero root of the diquark-antidiquark potential leads 
to smaller tetraquark masses which are shown with bold 
numbers in Table HE even with the same regularization 
cutoff the reported masses in Ref. [H are different from 
our results. 
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We should mention that in our calculations we have 
used 100 mesh points of Gauss-Legendre quadrature for 
integration over angle between r and r' with a linear 
mapping, whereas the integration over r' is done by a 
hyperbolic-linear mapping with the sub-intervals [0, -I- 

+ [rc, 20] GeV“^ using 150 mesh points to achieve 
the converged masses with four significant digits . 


TABLE II. Masses of charm diquark-antidiquark for Is state 
with total angular momentum J = 0 as a function of regular¬ 
ization cutoff Tc- The bold numbers are calculated for values 
of Tc which are non zero root of the diquark-antidiquark po¬ 
tential and numbers in parentheses are corresponding masses 
from Ref. [![. 


rc [GeV-^l 


Mass [GeV] 




cqcq 



SS 


AA 

2.831 

3.858 


3.987 (3.83908) 

2.844 

3.859 (3.70314) 


3.987 

2.85 

3.859 


3.988 

2.9 

3.862 


3.991 

3.0 

3.867 


3.995 

3.5 

3.879 


4.007 

4.0 

3.883 


4.011 

5.0 

3.885 


4.013 

6.0 

3.885 


4.013 



cscs 



SS 


AA 

2.808 

4.092 


4.226 (4.09962) 

2.821 

4.093 (4.05390) 


4.227 

2.85 

4.094 


4.228 

2.9 

4.097 


4.231 

3.0 

4.101 


4.235 

3.5 

4.112 


4.245 

4.0 

4.116 


4.248 

5.0 

4.117 


4.250 

6.0 

4.117 


4.250 


Although the conclusions based on the incorrect re¬ 


sults are misleading but we don’t understand this part of 
their conclusion: ” our results are in good agreement with 
the results derived from complicated relativistic methods 
and can be a good replacement for them”. How a nonrel- 
ativistic spin-independent method can be an alternative 
to a relativistic spin-dependent method, even when the 
results of both methods be in good agreement? 

In conclusion, our calculations demonstrate that all of 
the results published in recent letter by M. Monemzadeh 
et al. [l|, for masses of tetraquarks with total angular 
momentum equal to zero, are incorrect and the letter has 
serious numerical problems. The correct results should 
be as results given in Table Hi Moreover, beside these 
incorrect results, we have difficulty to understand how 
the authors have obtained the masses of tetraquarks for 
non zero total angular momentum in a spin-independent 
formalism! 

ACKNOWLEDGMENTS 

The author would like to thank V.O. Galkin for pro¬ 
viding the values of strong coupling constant for different 
tetraquark states. This work is supported by National 
Science Foundation under contract NSF-PHY-1005587 
with Ohio University. Partial support was also provided 
by the Institute of Nuclear and Particle Physics at Ohio 
University. 


[1] M. Monemzadeh, N. Tazimi, P. Sadeghi, Phys. Lett. B 
741, 124 (2015). 

[2] R. A. Malfiiet and J. A. Tjon, Nucl. Phys. A 127, 161 
(1969). 

[3] D. Ebert, R. N. Faustov, V. O. Galkin, Phys. Lett. B 634, 
214 (2006). 

[4] D. Ebert, R. Faustov, V. Galkin, Phys. Lett. B 696, 241 

( 2011 ). 

[5] M. R. Hadizadeh, A. Khaledi-Nasab, in preparation. 









